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ABSTRACT 

(N . We compare three different models of rotating neutron star spacetimes: i) the 

Hartle-Thorne (1968) slow-rotation approximation, keeping terms up to second order 
in the stellar angular velocity; ii) the exact analytic vacuum solution of Manko et al. 
(2000a, 2000b); and iii) a numerical solution of the full Einstein equations. In the 
first part of the paper we estimate the limits of validity of the slow-rotation expan- 
■ sion by computing relative errors in the spacetime's quadrupole moment Q and in the 

T^lj- I corotating and counterrotating radii of Innermost Stable Circular Orbits (ISCOs) R± . 

. We integrate the Hartle-Thorne structure equations for five representative equations 

of state. Then we match these models to numerical solutions of the Einstein equa- 
tions, imposing the condition that the gravitational mass and angular momentum of 
I i the models be the same. We find that the Hartle-Thorne approximation gives very 

^jr.' good predictions for the ISCO radii, with R± accurate to better than 1% even for 

the fastest millisecond pulsars. At these rotational rates the accuracy on Q is ~ 20%, 
and better for longer periods. In the second part of the paper we focus on the exte- 
rior vacuum spacetimes, comparing the Hartle-Thorne approximation and the Manko 
analytic solution to the numerical models using Newman-Penrose (NP) coordinate- 
' independent quantities. For all three spacetimes we introduce a physically motivated 

'quasi-Kinnersley' NP frame. In this frame we evaluate a quantity, the speciality in- 
dex S, measuring the deviation of each stellar model from Petrov Type D. Deviations 
from speciality on the equatorial plane are smaller than 5% at star radii for the faster 
rotating models, and rapidly decrease for slower rotation rates and with distance. We 
find that, at leading order, the deviation from Type D is proportional to (Q — Qkctt)- 
Our main conclusion is that the Hartle-Thorne approximation is very reliable for most 
astrophysical applications. 
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1 INTRODUCTION 

Pulsars are observed to rotate with periods as small as 1.56 
ms (Backer et al. 1982, Kulkarni 1995, Chakrabarty et al. 
2003). Fast rotation is also very important in newly-born 
neutron stars, which may undergo secular and dynamical in- 
stabilities (Andersson & Kokkotas 2001, Kokkotas & Ruoff 
2002). An accurate modelling of the spacetime of rotating 
stars is therefore important both in astrophysical studies 
and in investigations of sources of gravitational waves. The 
aim of this paper is to compare two relativistic approximate 
models, the slow rotation expansion of Hartle & Thorne 
(1968, henceforth HT) and the exact vacuum solution by 



Manko et al. (Manko 2000a, 2000b), with a fuU general- 
relativistic numerical one, quantifying the range of validity 
of the HT expansion and studying the possible overlap of 
the two approximations. 

The Schwarzschild metric describes either a non- 
rotating black hole (BH) or the exterior of a non-rotating 
star: this is a well-known consequence of Birkhoff's theo- 
rem. In contrast, the rotating Kerr BH metric describes a 
rotating star only to linear order in rotational velocity. At 
higher orders, the multipole moments of the gravitational 
field created by a rapidly rotating compact star differ from 
those of a BH. This difference in the multipolar structure 
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has important consequences for the observation of the elec- 
tromagnetic and gravitational radiation from these objects. 
For example, the detection of gravitational waves emitted 
by particles in orbit around BHs and compact stars can be 
used, in principle, to map the multipolar structure of the 
corresponding spacetime and check the validity of the BH 
no-hair theorem (Ryan 1995, 1997; Collins & Hughes 2004). 

A relatively simple description of the spacetime of ro- 
tating stars was developed many years ago by Hartle and 
collaborators (Hartle 1967, Hartle & Thorne 1968; see also 
Chandrasekhar & Miller 1974) using a slow rotation expan- 
sion, that is, an expansion in the parameter e = Here 
Q, is the star's angular velocity and Q,* = {M / is a 
rotational scale set by the Keplerian frequency of a test par- 
ticle sitting on the equator of a non-rotating star with grav- 
itational mass M and radius R (throughout this paper we 
use geometrical units, G = c = 1). In the Newtonian limit, 
uniformly rotating polytropic stars shed matter at the equa- 
tor at a mass-shedding frequency which is roughly given by 

^ (2/3)^/^ r2*. This is not true in a general relativis- 
tic framework (see eg. Stergioulas 2003). However Q,* does 
give an order of magnitude estimate of the mass shedding 
frequency, so we can roughly expect the HT treatment to 
break down when e ~ 1. 

The HT method allows a systematic computation of the 
interior and exterior structure at successive orders in e. At 
order e, deviations from non-rotating stellar models are due 
to the relativistic frame dragging, but the stellar structure 
is still spherical. Deviations from sphericity are measured by 
the star's quadrupole moment; they appear at order e^. Be- 
ing a slow-rotation expansion, the HT metric is expected to 
lose accuracy for rapid spin rates, but the limits of validity of 
the approximation are not well understood. Weber & Glen- 
denning (1992) compared the slow- rotation expansion with 
early numerical results by Friedman et al. (1986). Their com- 
parison led to 'compatible results' down to rotational Kepler 
periods Pk — 0.5 ms. However their self-consistency equa- 
tion [Eq. (7) in Weber & Glendenning (1992)] to determine 
the mass-shedding frequency is not consistently truncated 
at second order in angular velocity, as it should be. This is 
presumably one of the reasons their method becomes very 
inaccurate close to the limiting-mass model. 

Numerical solutions for uniformly rotating stars can 
now be obtained in both the slow and fast-rotation regimes. 
The HT approximation is quite simple to implement: our 
own numerical code will be described in Section |^ In con- 
trast, a solution of the Einstein equations for rapidly rotat- 
ing models is non-trivial. Komatsu, Eriguchi and Hachisu 
(1989) devised a self-consistent method to solve the full Ein- 
stein equations, which was later improved by Cook, Shapiro 
and Teukolsky (1994, henceforth GST). Stergioulas' RNS 
code (Stergioulas & Friedman 1995) is a variant of the GST 
code which has been extended to allow for strange-star mat- 
ter equations of state (EOSs), differential rotation, computa- 
tion of multipole moments and so on. A very accurate spec- 
tral method for computing rapidly rotating stellar models 
was independently developed by Bonazzola et al. (cf. Bonaz- 
zola et al. 1993, Salgado et al. 1994) and became part of 
the numerical relativity library LORENE (Bonazzola et al. 
1998). Nozawa et al. (1998) showed that the different nu- 
merical methods are in remarkable agreement, typical dif- 
ferences being of order 10~^ or smaller for smooth EOSs. 



A new multi-domain spectral method has been introduced 
by Ansorg et al. (2002), who achieved near-machine accu- 
racy for uniform density stars; their results are in very good 
agreement with those of the other groups. The RNS code 
therefore provides a numerically accurate description of fast- 
rotating relativistic stars. 

Berti & Stergioulas (2003, henceforth BS) used the RNS 
code to compute constant rest-mass (so-called evolutionary) 
sequences of three types using five realistic EOSs. The first 
type mimics the evolutionary track of a canonical neutron 
star having gravitational mass M = 1.4M0 in the non- 
rotating limit. The second type has the maximum grav- 
itational mass M = Mmax allowed by the given EOS in 
the non-rotating limit. In the following, we will sometimes 
refer to these sequences as type '14' and type 'MM'. The 
third type is supramassive: its baryon mass is larger than 
allowed by the non-rotating Tolman-Oppenheimer-Volkoff 
(TOV) equations. A supramassive sequence does not termi- 
nate with a non-rotating model, so that a comparison with 
HT models is not always possible. Therefore we do not con- 
sider supramassive models in the following. 

Another method of approximating the exterior space- 
time takes the form of a multipole expansion far from the 
star. This technique has been used by Shibata & Sasaki 
(1998) to derive approximate analytic expressions for the 
location of the Innermost Stable Circular Orbits (ISGOs) 
around rapidly rotating relativistic stars. However the tech- 
nique cannot easily be applied to other problems. Further- 
more, the relevant multipole moments must be computed 
numerically for each model one wants to analyze, and then 
plugged into the relevant formulas. 

Approximate analytic models for fast rotating stars 
which include only the most relevant multipolar structure 
of the exterior spacetime (that is, the lowest-order mass and 
current multipoles) would be very useful for astrophysical 
applications. In this paper we focus on an interesting fam- 
ily of solutions found by Manko et al. (2000a, 2000b). They 
used Sibgatullin's method (SibgatuUin 1991), which has the 
advantage that the multipolar structure of the spacetime 
can be prescribed a priori, to find an exact, nine-parameter, 
analytic solution of the Einstein-Maxwell equations. Consid- 
eration of the full Einstein-Maxwell system is motivated by 
the fact that compact stars normally have extremely large 
magnetic fields, which could significantly change their prop- 
erties (Bocquet et al. 1995) and the development of instabil- 
ities (Spruit et al. 1999; RezzoUa et al. 2000). However, for 
simplicity, we do not consider effects due to magnetic fields 
in this paper. 

An approximate description of a rapidly rotating neu- 
tron star can be obtained using a five parameter sub-case of 
the analytic solution (Manko 2000b). Ignoring charge and 
magnetic moment the free parameters reduce to three. This 
three-parameter solution is an exact solution of the Einstein 
equations in vacuum, but it only provides an approximate 
description of a neutron star exterior'^. 

^ BS focused on four out of the a priori infinite multipole mo- 
ments of the Manko metric: gravitational mass M, angular mo- 
mentum J, mass quadrupole moment Q and current octupole 
moment 53. However their statement that these are the only 
nonvanishing moments is not precise: there do exist higher-order 
nonvanishing physical moments of the Manko metric. 
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BS showed that the Manko solution approximates 
rather well the field of a rapidly rotating neutron star in 
the fast rotation regime. They used the three parameters of 
the analytic solution to match the dominant multipole mo- 
ments (M, J and A/2 = —Q) to given numerical models. 
They found that an indicative measure of the deviation of 
the Manko metric from the numerical metric is given by the 
relative error in the first non-freely adjustable multipole mo- 
ment: the current octupole S3. BS found a peculiar feature of 
the Manko solution: the solution's quadrupole moment can- 
not be matched to the quadrupole moment of numerically 
computed spacetimes for slow rotation rates. Therefore, it 
is of interest to understand whether the regimes where the 
HT and the Manko metrics are applicable are completely 
distinct, or have some range of in common. 

This paper compares the HT slow rotation expansion 
and the exact vacuum solution by Manko et al. with the 
full general-relativistic numerical CST spacetime. We will 
assume that our numerical models, obtained using the RNS 
code, are exact (within their numerical accuracy, that is of 
the order of 10"'' in each computed quantity). An overview 
of the three different spacetimes we consider in this paper is 
presented in Section|51 where we also introduce our notation. 

The rest of the paper consists of two parts. In the first 
we integrate the HT equations of structure for the same set 
of 'realistic', nuclear physics motivated EOSs examined in 
BS (Section From the stable models in the HT sequence 
we choose those which have the same gravitational mass 
and angular momentum as a given numerical model. The 
procedure to do this is developed in Section [3. II 

In Section^we quantify deviations of the HT and CST 
spacetimes. A trustworthy measure of the deviation is given 
by the relative error in the first non-freely adjustable mul- 
tipole moment in the HT spacetime: M2 = — Q. Picking a 
representative EOS we show that deviations in Q are at most 
~ 20% at the rotation rates corresponding to the fastest ob- 
served millisecond pulsar PSR J1939-I-2134 (Backer et al. 
1982), for both mass sequences. We also use some recent 
results by Abramowicz et al. (2003) to compute the ISCO 
radii predicted by the HT approximation for sequences of 
slowly-rotating models. We find an interesting result: the 
HT approximation gives an excellent approximation to the 
ISCO predicted by numerical solutions, relative errors being 
smaller than 1 % for astrophysically relevant rotation rates, 
down to the 1.5 ms period of PSR J1939-I-2134. 

In the second part of the paper we compare the three 
vacuum exterior spacetimes, focusing on invariant quanti- 
ties, and looking at their Petrov classification. In a nutshell, 
the Petrov classification of a spacetime depends on which of 
the Weyl scalars can be set to zero by suitable tetrad rota- 
tions (Chandrasekhar 1983, Stephani et al. 2003, Stephani 
2004). The Kerr metric is of Petrov Type D, hence alge- 
braically special. In general, rotating perfect fiuid spacetimes 
will not be of Type D (see eg. Fodor and Perjes 2000), but 
will be of (general) Type I. Our invariant comparison of the 
three different models involves quantifying the deviation of 
each metric from speciality. Indeed one might wish to as- 
sume (for example in perturbative calculations of gravita- 
tional wave emission) that rotating neutron star spacetimes 
are in some sense 'close' to Type D. Then one should de- 
velop a tool to quantify their 'non-Type-D-ness': here we 
use the curvature invariant speciality index S introduced by 



Baker & Campanelli (2000) . Each metric is expressed in dif- 
ferent coordinates, so we evaluate S at the stellar radius, 
a physically identifiable radial location. This quantity may 
also indicate whether it is sensible to carry out perturba- 
tive calculations on neutron star backgrounds in a 'Type-D- 
approximation', and to extract the gravitational wave con- 
tent of the spacetime through a Teukolsky-like formalism. 
In this paper, however, we will focus only on calculations of 
the Weyl scalars in the background spacetimes. 

We choose the simplest and most physically transpar- 
ent tetrad to compute the Weyl scalars and to characterise 
the non-Type-D-ness of the spacetime. First of all we con- 
struct a transverse tetrad, in which = i/)3 = 0, so getting 
rid of this tetrad gauge freedom (Szekeres 1965). Then we 
carry out a tetrad rotation to go to the canonical frame in 
which ipo = ip4 = V'04 (PoUney et al. 2000; Re et al. 2003), 
that we call symmetric. In a Type I spacetime, there are 
three distinct equivalence classes of transverse frames (Bee- 
tle & Burko 2002) , each containing a symmetric tetrad. From 
these three different symmetric tetrads we select a 'quasi- 
Kinnersley' tetrad (cf. Nerozzi et al. 2004), i.e. the one that 
directly reduces to the canonical frame for Type D when we 
set the parameters to reduce the metric to this type. 

We apply this procedure, which is described in detail 
in Section |S] to all three metrics. In Section 15] we show the 
scalars and the corresponding speciality index. We find that 
deviations from speciality on the equatorial plane are smaller 
than 5% at the stellar radius, even for the faster rotating 
models. These deviations rapidly decrease for slower rotation 
rates and at larger distances from the star. We find that, at 
leading order, the deviation from Type D is proportional to 
{Q — Qkctt). The conclusions and a discussion follow. 



2 THE THREE METRICS 

The three metrics we will consider in this paper (HT, Manko 
and CST) can be written in the following general axisym- 
metric stationary form, with coordinates [t, ci, C2, (j}] 

ds^ = goodt^ + 2go3dtd(t> + gsadcj}'^ + giidcl + g22dcl (1) 

where all the metric coefficients are functions of ci and C2 
only. Being stationary and axisymmetric the spacetime has 
two Killing vectors, one timelike — [1,0,0,0] and one 
spacelike ip" = [0,0,0,1]. 

The metrics can be explicitly written as: 

(ds^)HT = -e''{l + 2h)dt^ + e^[l + 2m/{r~2M)]dr^ 

+ r^{l + 2k)[de^ + sin e[d(j> - ujdtf] (2) 

{ds^)Manko = f^^[e^'^{dp^ + dz^)+p^d(j)^] 

- f{dt-Cbdij>f (3) 

{ds'^)csT ~ —e^^dt? + siri^ 6{d(j) — Ludt)^ 

+ e^"{dr^ + f^de'') (4) 

Ignoring terms of 0{e^), the HT metric components are 
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(cf. Abramowicz et al. 2003)^ 

9rr = (^1 - ^) [1 + fGi + + O(e^) 
g,, = - [1 - [1 + /Fi - qF2] + 0{e') 

gee = r 
g<t>4> = 
9t,j, = 

where 



, sin e + 0{e) 



sin^6l + 0(e^) 



Fi = -pW + Ai, 

F2 = 5r^p{3u^ - 1) - M) {2Af + 6Mr - 3r^) - Ai, 



Al = 



15r(r-2M)(l-3M^) 
16M2 



In 



r - 2M / ' 

r 



15(r^ - 2M^)(3m^ - 1) 
^ ~ 16M2 ^\r-2M^ 

Gi = p[(L-72MV)-3«^(L-56MV)] -Al, 

Hi = A2 + {8Mr*)-\l-3u'^) X 

X (16M^+8M'*r-10M^r-^ + 15Mr'* + 15r^), 

H2 = -A2 + (8Af-r)-^5(l-3u^)(2M^-3Mr-3r2), 

L = (80M*' + 8Af "^r^ + IOMV^ + 20M V - 45Mr^ + 15r® 

p = [8M/(r-2M))]"\ 

W = (r-M)(16M^+8M*r-10MV^-30Mr* + 15/') 

+ w^(48M®-8MV-24M*r^-30M^r^ 

- 60MV* + 135Mr^ - 45r^), 

with u = COS0. The dimensiorilcss angular momentum and 
quadrupole moment j and q are defined as: 



j = J/M^ 



q = Q/M^ 



(5) 



The fill! Manko metric is rather lengthy and can be 
found in many papers (eg. Manko et al. 2000a, 2000b; Stute 
& Camenzind 2002; BS), so we will not give it explicitly 
here. We only recall that it depends on three parameters: 
the gravitational mass M, the angular momentum per unit 
mass a = J/M and a third parameter b which is related to 
the quadrupole moment by 



Q = -M{d-S-ab + a^), 



(6) 



where 
5 = 

d = 



1 

4 [M2 - (o - 6)2] ■ 



The next non-zero multipole moment for the Manko metric 
is the current octupole, which is given by 



53 = -M [a^ - 2c?h + o[6^ + 2(d - S)\ - h(d - <5)] 



(7) 



2 The original Hartlc (1967) metric differs from tfie one presented 
here, in that it was not consistently truncated to order e2. The 
Abramowicz et al. (2003) preprint contains some minor sign er- 
rors. 



For a detailed description of the CST metric we refer 
to the original paper (Cook et al. 1994). The extraction of 
multipole moments from the CST metric and the numeri- 
cal calculation of ISCOs along evolutionary sequences are 
discussed in BS. 

In vacuum, the HT metric is equivalent to the Kerr 
metric at second order in J when one sets Q = /M, and 
reduces to the Schwarzschild solution when J = Q = 0. 
It is not possible to obtain Kerr directly from the Manko 
metric by a smooth variation of the real parameters a, M 
and 6. In other words, Kerr and Schwarzschild are not part 
of the 3 real parameter family of solutions described by the 
Manko metric. However they can be obtained by analytic 
continuation, setting = a^—M^. Schwarzschild is a special 
case with a = and 6^ = —M^. 



3 INTEGRATION OF THE HARTLE-THORNE 
EQUATIONS OF STRUCTURE 

In this Section we outline some important aspects of the HT 
equations of structure. The relevant formalism was originally 
developed in (Hartle 1967), and the integration method is 
summarized in (Hartle & Thorne 1968). In the numerical 
code we used the reformulation of the equations of structure 
given in Appendix A of (Sumiyoshi et al. 1999). We refer to 
that paper for the equations ^. 

' As discussed in the Introduction, the HT metric is a 
slow rotation expansion in the angular velocity truncated to 
second-order. We consider the expansion parameter to be 



(8) 



The solutions of the TOV equations for a given EOS form 
a one-parameter family, where, for example, one can use 
the star's central energy density pc as the parameter. Simi- 
larly, solutions of the HT equations of structure form a two- 
parameter family. Natural parameters to label each model 
are the central energy density pc and the rotational param- 
eter e. 

Arriving at a particular HT model involves three con- 
ceptual stages. First, the stellar structure of a non-rotating 
star is obtained from the TOV equations, giving the non- 
rotating gravitational mass M and stellar radius R. Note 
that these non-rotating quantities are unlabellcd. From 
these the rotational scale O* = ^ M/R^ can be calculated. 
We can also compute the star's binding energy Eb, defined 



^ Comparing with the HT original papers, we found that Ap- 
pendix A of (Sumiyoshi et al. 1999) contains a few typos: 

• In Equation (All), r** should be replaced with r^^. 

• In Equation (A14), i?"* should be replaced with r*. 

• In the last factor of Equation (A15), r{r — 2m) should be 
replaced with (r — 2m). 

• In Equation (A36), p should be replaced with pc. 

• In Equation (A41), e~''/2 should be replaced with e"/'^. 

• In Equation (A54), go should be replaced with i;2. 

Notice also that Sumiyoshi et al. do not explicitly give boundary 
conditions at the center for the 'particular' solutions fej: ^2 ■ These 
boundary conditions are given in equations (128)-(130) of Hartle 
(1967); however, in equation (130) of Hartle's paper the fa<;tor 4it 
should be replaced by 2it. 
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as the difference between the baryonic mass Alt and the 
gravitational mass M: 



Eb = Mb 



M. 



(9) 



Second, a model with e = 1, i.e. a model rotating at the 
'maximum' angular velocity of fi*, is obtained from the 
full HT equations. Here all quantities are labelled with 
a superscript asterisk : the model has gravitational mass 
M* = M + 5M* , angular momentum J*, and quadrupole 
moment Q* . Finally, the model with the chosen value of 
e < 1 is calculated, having rotating gravitational mass 
M''^' — M+SM^'^\ angular momentum and quadrupole 
moment Q^'^K 

The effects of rotation on the various quantities of in- 
terest appear at different orders in e. The general relativistic 
frame dragging appears at first order, and the angular mo- 
mentum is; 



(10) 



The radius of gyration Rg, defined so that the moment 

2 



of inertia is 7 = 
J \ 1/2 



MRg, can then be calculated as 



Rn 



M 



\nM ) 



1/2 



J* 

TFm 



1/2 



+ (11) 



where the first two relations are exact, and the third ex- 
presses the fact that computing Rg requires the integration 
of the equations of structure at first order in e. However, 
involving the ratio of J and Q. (which are both proportional 
to e) Rg does not depend on e. For the same reason, com- 
puting the second order contribution to Rg would require 
third order corrections to J and Q.. 

At second order in e there appear corrections to the stel- 
lar eccentricity Cs, quadrupole moment Q, mass M, binding 
energy Eb, and coordinate stellar radius R. The stellar ec- 
centricity is 



= [{Re/Rp) 



,1/2 



(12) 



where Re and Rp are the stellar coordinate radii at the equa- 
tor and at the pole, respectively*. The quadrupole moment 
Q can be invariantly defined as the leading coefficient in the 
asymptotic behavior at large distances from the origin of 
the eccentricity of a family of spheroidal surfaces where the 
norm of the timelike Killing vector is constant. Taking into 
account that the zero order (non-rotating) values of Ca and 
Q are zero, we have 



(13) 



This means, in particular, that for a fixed value of the 

* In the notation of Hartle & Thorne (1968) the stellar radius 
R(e) = R + gQ(ij) + ^2(i?)P2(cos6'), where R is the non-rotating 
stellar radius, ^o{'>') and ^2{t) are radial functions determined 
by the HT perturbation equations and P2 (cos 6) is a Legendre 
polynomial. The equatorial radius Re = R{tt/2), and the polar 
radius Rp = R{0). 





2 * 












= e^SE-B 


sr'^'^ 


= e^SR*. 



central energy density, a stellar model rotating with angular 
velocity O = eSl* has a gravitational mass 



M 



(14) 



We have checked our numerical code by reproducing the 
results shown in Hartle & Friedman (1975) for polytropic 
EOSs. We found agreement on all significant digits for the 
quantities we have just defined. Then we chose the same 
set of 'realistic' EOSs used in BS, and verified that in the 
non-rotating limit our numerical results are consistent with 
Stergioulas' RNS code (described in Stergioulas & Friedman 
1995). 

For any given 'realistic' EOS, we can construct a se- 
quence of rotating models by integrating the HT structure 
equations for a series of values of the central energy den- 
sity. Our results are summarized in Tables [TIKI The chosen 
central energy density pc, radius R, gravitational mass M, 
binding energy Eb and fi* are listed in the first five columns 
(marked by e", since they do not depend on rotational cor- 
rections). At order e we list the value of the radius of gy- 
ration R*g (from which the angular momentum J* can be 
easily obtained) and of the frame dragging function uie at 
the center (normalized by fi*). Finally we list the second- 
order quantities 6R*/R, SM*/M, SEb/Eb, Q* /{MR^) and 
e*. 

Let us stress again that although f2* is not the mass- 
shedding limit in a general relativistic framework, because of 
its physical meaning it does give an order of magnitude esti- 
mate of the mass shedding frequency, and so we can expect 
the HT treatment to break down when e ~ 1. 



3.1 Matching the Hartle-Thorne and numerical 
spacetimes 

To find a HT model that corresponds to a particular GST 
stellar model, say one belonging to the constant-rest mass 
sequences computed in BS using the RNS code, we adopt 
the following procedure, in which quantities obtained from 
the GST model are labeled with overbars : M, J and Q. 

The aim is to find a HT model for which A/'"' = M and 
J^"^' = J. Imposing these values on Eqs. HUH and 11411 gives 



J/J* 



and 

M = 



M + {J/rf5M\ 



(15) 



(16) 



In 1161 1. the quantities J*, M and 5M* are effectively func- 
tions of pc, since once pc is specified they can be calculated 
from e = or e = 1 models. Therefore, in order to find the 
matching HT model, we must search for a value of pc for 
which Eq. 1161 holds; e is then given by Eq. 1151 . 
In practice, we consider the function 



f[pc) = M- [M + [J/rfSM*] 



(17) 



To locate the matching model with a given accuracy (that 
we set to 10~^) we search for the root of Eq. 1171 using 
bisection. Note that, since the HT expansion is a small-e 
approximation, we expect to find solutions only for small 
values of the angular momentum J. This turns out to be 
indeed the case. 

In Table|S|we show some basic features of the matching 
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models for the five EOSs we consider. The HT models we 
compute are matched to the slowest-rotating models of the 
evolutionary sequences presented in Tables 1-5 of BS: as just 
mentioned, Eq. (I17II has no solutions beyond some critical 
rotation rate. 

For each EOS we present matching models for the 14 se- 
quence above matching models for the MM sequence. From 
left to right we list: the central energy density pc for which 
we find a matching solution (if such a central energy density 
exists); the value of e = J / J* at the given central energy 
density; the gravitational mass M and the angular momen- 
tum J, in geometrical units; the value of the quadrupole 
moment predicted by the HT equations of structure for the 
selected model, which scales as e^, and is given by 

-M2 = Q<'^ = (J/J*)'Q*; (18) 

and finally, 

5Q = 100(O''> - Q)/Q, (19) 

that is, the percentage deviation of the HT quadrupole mo- 
ment from its value Q for the numerical spacetime, which 
can be found in Tables 1-5 of BS. 



4 HOW ACCURATE IS THE 

SLOW-ROTATION EXPANSION? 

In the previous Section we have seen how, given a numer- 
ical solution of the full Einstein equations with gravita- 
tional mass M = M and (small enough) angular momen- 
tum J = J, we can construct a HT model matching the 
given values of M and J. Our purpose in this Section is to 
give two alternative, quantitative estimates of the difference 
between each of the HT matching models and the corre- 
sponding numerical model. The first estimate will be based 
on the calculation of the spacetimes' quadrupole moment as 
a function of the rotation rate; the second, on the calculation 
of the ISCO radii. 

4.1 Deviations in the quadrupole moment 

The matching models described in Section 13.11 are con- 
structed by imposing that the first mass multipole Af and 
the first current multipole J be the same as in the numeri- 
cal spacetime. The HT solution is a two-parameter family: 
once pc and e are specified, the whole multipolar structure 
of the HT spacetime is. In particular, at second order in 
e the HT structure equations predict a specific value for 
the quadrupole moment Q^^. Since multipole moments are 
global features of a spacetime (Fodor, Hoenselaers & Perjes 
1989), we can measure the 'distance' of the HT spacetime 
from the numerical spacetime by looking at the difference 
between the value Q^^ predicted by the HT code, and the 
'exact' value Q predicted by the RNS code. That is, we com- 
pute the relative deviation 

AQ _ Q"^~Q 

Q ' Q ■ ^ ' 

Given that the HT metric is exact in the limit e ^ 0, but 
becomes less accurate for larger values of e, we expect this 
deviation to increase with the rotation rate: an explicit cal- 
culation shows that this is indeed the case. 




0,0 I ■ ' ■ ' ■ ' ■ ' ■ <— 

0,0 0,1 0,2 0,3 0,4 0,5 

s 



Figure 1. Relative error in the quadrupole moment for the FPS 
14 sequence (squares) and the FPS MM sequence (circles). 

As shown in Table |S| matching models - that is, so- 
lutions of Eq. (1171 - only exist up to some critical value 
of the angular momentum J (and of the expansion param- 
eter e) where the HT approximation breaks down, and is 
unable to reproduce the numerical results. The deviation of 
from the 'true' numerical value when e ~ 0.2 depends 
slightly on the EOS, but it is typically ~ 10% for the 14 
sequence and ~ 20% for the MM sequence. For the 14 se- 
quence, the relative error becomes larger than 20 % when 
e ~ 0.4. From Table |S1 one can see that deviations are not 
very sensitive to the EOS: for example, when e ~ 0.4 the rel- 
ative error for the MM sequence ranges (roughly) between 
25 % and 30 %. At fixed e, AQ/Q is larger for the MM se- 
quence than for the 14 sequence. However recall that Q is 
generally larger for low-mass evolutionary sequences, where 
gravity is weaker and centrifugal effects give the star a more 
oblate shape. Deviations in the ISCO radii, as given by Eq. 
1211 1 below, depend on the absolute value of Q. Since Q is 
larger for the 14 sequence, and percentage errors induced 
by the HT approximation are comparable to those for the 
MM sequence, the HT approximation should induce larger 
errors in the ISCO for the low-mass 14 sequence. This will 
be explicitly shown in Section f4.2l 

Table 6 contains only a few slow-rotating models for 
each EOS. To get a better idea of what happens in the slow- 
rotation regime we have selected EOS FPS (which is quite 
average in terms of stiffness) and computed more numeri- 
cal models for small values of J (that is, for small values 
of e). Namely, we used the RNS code to compute 11 con- 
stant rest-mass models for the 14 sequence and 10 models 
for the MM sequence, spanning the range e = to e ~ 0.5. 
To each of these models we matched a HT model. Results 
for the relative deviations in the quadrupole moment are 
shown in Fig. The convergence to zero as e — » is very 
smooth for the FPS 14 sequence. For the FPS MM sequence, 
the slowest rotating model seems to show a slight deviation 
from the general trend. We suspect this is due to a loss of 
accuracy of the RNS code in computing the comparatively 
smaller Q (for any given e) of the more compact MM star 
models with small rotation rates. Indeed, when we compute 
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the curvature invariants for these models we also find some 
numerical noise: this can be seen by a close scrutiny of the 
bottom curve in the lower panels of Figs. |^ |^ and The 
relative error AQ/Q is of the same order of magnitude in 
both cases, but it is systematically larger for the FPS MM 
sequence than for the FPS 14 sequence. 



4.2 Deviations in Innermost Stable Circular 
Orbits 

The exterior vacuum metric for a HT model truncated at 
second order in e depends only on M, J and Q = Q^^ . In 
particular, once we have computed these multipole moments 
we can use them to compute the coordinate ISCO radius 
using the results by Abramowicz et al. (2003). Using the 
dimensionless quantities j and q introduced in Eq. @, their 
prediction for the coordinate ISCO radii is 



IT. (3 



3/2 



.2 /251647 3 



9325 3 
-^ + 2401n- 



(21) 



where r+ refers to corotating orbits, and r_ to counterrotat- 
ing orbits. The comparison of the HT ISCO radii as given 
by Eq. I)21|l to numerical results from the RNS code provides 
an alternative measure of the difference between a numerical 
model and the corresponding HT matching model. 

We consider again the slow-rotating evolutionary se- 
quence for EOS FPS described in Section 14.11 For each 
model along this evolutionary sequence we compute M, J 
and Q^^, and we plug those values in Eq. (I21t to compute 
the coordinate ISCO radii. Then we can easily compute the 
circumferential ISCO radii for these equatorial orbits as 



r±,e = n/2), 



(22) 



and compare our results with the numerical value of the 
circumferential ISCO radii^ as obtained from RNS. 

In Fig.|21we show results for the FPS 14 sequence (top 
panel) and for the FPS MM sequence (bottom panel). In 
each panel we plot together, as functions of e: 1) the cir- 
cumferential radii R of the stars, 2) the corotating ISCOs 
R+, 3) the counterrotating ISCOs R-. Filled symbols are 
obtained from RNS, and empty symbols correspond to the 
HT matching models. The circumferential radius for each 
matching model is just 

(HT) 



R 



\/ g4.4,{Re,0 = 7r/2), 



(23) 



where Re is the coordinate equatorial radius, obtained by 
integrating the HT equations of structure. 

It is evident from Fig.|5]that for neutron stars of canon- 
ical mass M ~ I.4M0 the deviation of HT from the full nu- 
merical solution is larger. However, even in that case the HT 
approximation gives very accurate predictions both for the 
stellar equatorial radii and for the ISCO radii up to e ~ 0.3. 
In particular, the fastest rotating model in the FPS 14 se- 
quence for which a corotating ISCO exists has e = 0.318 
and a spin frequency in physical units J^phys = f2phys/27r = 



^ These are given by R± 
BS. 



Re + h± in the notation of CST and 



R. R,<,,Jkm) 




R, R,sco(km) 
20 




Figure 2. Filled circles give the stellar circumferential radius 
_R, filled diamonds the circumferential radius of corotating ISCOs 
_R+ (when it exists) and filled squares the circumferential radius 
of counterrotating ISCOs i?_ for the numerical spacetime. Empty 
symbols give the corresponding quantities for the HT metric. The 
top panel refers to the FPS 14 evolutionary sequence, the bottom 
panel to the FPS MM sequence. 



599 s~^, corresponding to an orbital period P = 1.67 ms. 
Therefore our calculation shows that the HT approximation 
accurately predicts ISCO radii down to periods which are 
comparable to that of the fastest observed millisecond pul- 
sar, PSR J1939-I-2134, spinning at P = 1.56 ms (Backer et 
al. 1982). 

Relative errors for corotating and counterrotating IS- 
COs are shown in Fig.Ejfor both the FPS 14 and FPS MM 
sequences. In this plot we use the same convention as in BS: 
e < corresponds to counterrotating orbits. As one would 
expect, errors tend to zero as e — > 0. The error for counter- 
rotating ISCOs monotonically increases for large e; the error 
for corotating ISCOs does not (at least for the FPS MM se- 
quence) and it is typically extremely small (below 5 parts 
in a thousand, roughly the same accuracy as our numeri- 
cal runs of the RNS code). We suspect that this could be 
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AR,sco/R,s 




AQ/Q 




Figure 3. Relative error in the ISCO radii, - R+)/R+ 

and (R^j"'^^ - R-)/R-, for the FPS 14 sequence (squares) and 
the FPS MM sequence (circles). In this plot we used the same 
convention as in BS: a negative e (and a negative J) corresponds 
to counterrotating orbits. 



due to the second order HT approximation being sufficiently 
accurate to compute tiie corotating ISCO, but not the coun- 
terrotating ones. To confirm our conjecture we would need 
to push the HT calculation to third order, which is beyond 
the scope of this paper. 

As anticipated in Section ISCO radii for the MM 
sequence agree better with the numerical results. Stars in 
the lower-mass FPS 14 evolutionary sequence have larger Q 
(they are more oblate), and relative errors in Q for both se- 
quences are roughly comparable. Therefore effects due to the 
quadrupole moment on the location of the ISCO are larger 
for the lower-mass sequence. Notice also that, although devi- 
ations in the quadrupole moment can be as large as 20 % for 
stars rotating with e ~ 0.3, the relative errors on the ISCO 
are typically smaller than 1 % at the same rotation rates. 
This can be understood by looking at Eq. 12H . For small ro- 
tation rates we have seen in Sectionj^that q — aj^ , where a 
is a constant of order unity. This scaling property is actually 
more general, applying also to numerical solutions with large 
rotation rates: Laarakkers & Poisson (1999) showed that a 
depends on the EOS and on the mass of the stellar model, 
ranging between a ~ 2 for high-mass models with soft EOSs 
and a ~ 12 for low-mass models with stiff EOSs (cf. their 
Table 7; see also Miller 1977 for similar results in the HT ap- 
proximation). The numerical coefficient of the j^-dependent 
correction in Eq. (I21II is ~ —0.226, while the numerical coef- 
ficient of the g-dependent correction is ~ 0.176: that is, the 
two coefficients are roughly equal in magnitude but oppo- 
site in sign. For this reason, shape-dependent effects in the 
ISCO (which are 0{e^)) tend to cancel for small rotation 
rates. This explains why the HT approximation does such a 
good job in predicting the ISCO radii. 

To better illustrate the physical significance of our re- 
sults, in Fig. 2] we plot the relative errors in the quadrupole 
moments (top panel) and in the ISCO radii (bottom panel) 
as a function of the spin period P (in milliseconds) of our 



^'^ISCO^'^ISCO 




Figure 4. The top panel shows the relative error in the 
quadrupole moment for the FPS 14 sequence (squares) and the 
FPS MM sequence (circles) as a function of the spin period P of 
the star (in milliseconds). The bottom panel shows the relative er- 
ror for corotating {R+) and counterrotating (R—) ISCO radii for 
the FPS 14 sequence (squares) and the FPS MM sequence (cir- 
cles) as a function of P. In both panels, the vertical line marks 
the spin frequency of the fastest observed millisecond pulsar PSR 
J1939-I-2134. 



rotating stellar models. This plot is useful for the follow- 
ing reasons. First of all, the relation between the HT ex- 
pansion parameter e and the rotational period of pulsars is 
non- linear. Not only is P proportional to 1/e, it also has (for 
the constant-rest mass sequences we consider) a weak but 
non-trivial dependence on the parameters M and R: 

2^^^^ 2^^ (24) 

fl eft" €{M/R^y/^ • ' 

Fig. |1| shows that all of our so-called 'slow rotation' mod- 
els are actually millisecond pulsars, many of which rotate 
even faster than the fastest observed millisecond pulsar PSR 
J1939-I-2134 (Backer et al. 1982). This gives a concrete idea 
of the rotation rates we are considering in this paper. 

In both panels we mark by a vertical line the spin fre- 
quency of PSR J1939-I-2134. The relative error in Q for a 
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pulsar spinning as fast as PSR J1939+2134 is of order 20 %, 
and this is quite independent of the stellar mass. Remark- 
ably, the relative error in corotating and counterrotating IS- 
COs induced by the HT approximation is smaller than 1 %, 
even when the pulsar rotates with a period as small as 1.5 
ms. This upper limit on deviations in the ISCO applies both 
to the FPS 14 sequence and to the FPS MM sequence. 

For concreteness, in Fig. 2] we presented the specific 
example of PSR J1939-I-2134. However we wish to stress 
that our analysis can be considered more general. Even if 
faster pulsars were discovered, our conclusions are likely to 
be only marginally altered. Indeed, Chakrabarti et al. (2003) 
performed a statistical analysis of observations of 11 nuclear 
powered pulsars with spin frequencies between 270 Hz and 
619 Hz. Their analysis implies an upper limit of 760 Hz (95% 
confidence level) on the spin frequency, to be compared with 
the 641 Hz of PSR J1939+2134. Therefore PSR J1939+2134 
should be very close to the maximum spin frequency allowed 
for an isolated neutron star. Interestingly, the upper limit 
from burst oscillations is much lower than the mass-shedding 
limit predicted by most EOSs: some mechanism must be at 
work to halt pulsar spin-up while accretion is still active. 
A plausible explanation is gravitational wave emission, due 
either to the excitation of the r-mode instability or to the 
presence of an accretion-induced crustal quadrupole moment 
(see Chakrabarti et al. 2003 for a discussion). 



5 WEYL SCALARS IN A 

QUASI-KINNERSLEY FRAME 

In the remainder of this paper we concentrate on the vacuum 
surrounding our rotating stellar models. We compute the 
Weyl scalars for the HT, Manko and CST spacetimes and 
define the 'quasi-Kinnersley' frame fSection l5.11 . In Section 
I5.2l we present the simple analytical results we obtained for 
the HT metric, and in Section lb.3l we describe the analogous 
calculation for the Manko and CST metrics. In Section|S|we 
give numerical results for all three spacetimes. 



5.1 The quasi-Kinnersley tetrad 

We introduce a Newman- Penrose null tetrad n", m", 
m"}, where and n° are real and m° and m" are complex 
conjugate, satisfying the usual orthonormality conditions: 
I'^Ua ~ 1, m'^fha = —1, all other products zero. The Weyl 
scalars are the usual contractions of the Weyl tensor. 



ipo = —Cpqrsl^rnTm^ , 

ipi = —Cpgrsl^nTm" , 

i>2 = —CpqrsV'rn'fn'n!' , 

tpi = —Cpqrsl^n'fnn 



(25) 



ypqr 

1p4 — —CpqrsTL^fn'n^Tfl' 



The speciality index (Baker & Campanelli 2000) is defined 



5 = 27JV/^ 



(26) 



where the curvature invariants I and J can be expressed as 



3(^2 



- 4l/)lt/j3 -I- ■lpo'4'4. 



We use the additional conventions that an NP tetrad with 
'/'I = ^3 = is called transverse (Beetle & Burko 2002), and 
that if further ijjo = tp4 then it is called symmetric, in which 
case we may refer to tpo = ip4 as tpoi- 

The interpretation of the Weyl scalars in the context 
of BH perturbation theory is that ipo represents radiation 
along the 1°" direction, 7/J4 represents radiation along the 
direction, while tp2 describes the Coulomb effect due to the 
presence of a central mass and the frame dragging of the 
background Kerr spacetime. ^1 and ips are longitudinal ef- 
fects of no physical interest that can be set to zero by the re- 
maining tetrad freedom (Teukolsky 1973, Stewart & Walker 
1974). In our case, however, and in general for a stationary 
spacetime of Petrov Type I, tpo and represent transverse 
curvature deviations (cf. Szekeres 1965) from the Coulom- 
bian ■!/;2 field due to rotation: see Eq. 1401 . 

It is useful then to consider only transverse frames. In 
a Type D spacetime a frame can be found in which ■ip2 is 
the only non-zero scalar: we refer to this as the Kinnersley 
frame (cf. Teukolsky 1973). With an appropriate choice of 
parameters (see below) both the Manko and the HT metrics 
can be reduced to Type D spacetimes (i.e. Schwarschild and 
Kerr). Here we work in the symmetric tetrad which, under 
these parameter choices, becomes the Kinnersley frame. In 
this tetrad tpo4 is an effect due solely to the non-Type-D-ness 
of the full metrics. We call this tetrad the 'quasi-Kinnersley' 
tetrad (Nerozzi et al. 2004). 

In a Type I spacetime, there are three distinct equiv- 
alence classes of transverse tetrads (Beetle & Burko 2002). 
Each class consists of all those transverse tetrads that can 
be related to each other by either the exchange freedom 
^ n", or by a class HI rotation®: 



I - 
n 



* An 



■02 



A- 
A- 

■02 



2„2i9 ; 



>1 



7/J4 



Ae-' 
A^e 



>3 
2i9„;, 



We find one particular transverse tetrad by making the 
ansatz that the real null vectors of the tetrad should be 
linear combinations of the Killing vectors ^" and 



= AC + Bv\ 



(27) 
(28) 



We then impose the NP tetrad inner product conditions and 
make particular (arbitrary) choices of the remaining free pa- 
rameters and signs (different choices may result in a different 
tetrad, but this would be related to the first by a class HI 
rotation, and so would be in the same equivalence class). 
The result is 



r = 



%/2 
1 . 



'533 



%/2 
1 
71 



1,0,0, 



1,0,0, 



(-ffo3 - 0-) 



g-Az 



933 



where 



= 903-900933- 



(29) 
(30) 
(31) 

(32) 



-{i'2)' + i'0t^'2lp4 + 2lpl1p2i'3 - 1p4{lpl) - i'oii'3) 



^ The parameter 9 used in this Section has obviously nothing to 
do with the angular coordinate in the HT metric. 
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To pick out a particular member of a class of tetrads one can 
choose the symmetric tetrad. This can always be found from 
a transverse tetrad through the following class III rotation 



'4>2 = ^2, 



A = |l/)o|/|V'4 



•^04 = l/'O = l/'4 = VV'0^4, (33) 



— — arctan 
4 



; (1/^0^4 



3? (V'oV'4 



(34) 



where the i/)'s are in the symmetric tetrad, the i/j's are in 
the non-symmetric tetrad and an overbar means complex 
conjugate. We do not write down the expressions for the 
tetrad vectors in the symmetric frame, since the parameters 
A and 6 are very long if written out in full, making the 
expression unenlightening, but it is worth noting that since 
this is a class III rotation, the directions of the and n" 
vectors do not change. 

We will call the above tetrad, built on the Killing vec- 
tors, Tl. Using a method suggested by Burko (2003) we can 
show that if the non-zero scalars in the symmetric version of 
Tl are ip^l and ip2^ then in the other two transverse sym- 
metric tetrads (T2 and T3) the non-zero scalars are given 
by: 



T2 



T3 



, T2 
V2 



, T2 



I Tl 



2 yVoA + ^V: 



§ ( -^04 



^ (^04 - 31/^ J^) 



(35) 
(36) 
(37) 
(38) 



Hence, once the scalars in one transverse tetrad are 
known, it is easy to calculate those in the other two frames. 
We have done this in the reductions of Manko and HT to the 
Type D metrics described in the next section. In each case 
it turns out that T3 is the Kinnersley tetrad, i.e. that frame 
in which, for the Type D metrics, il>2 is the only non-zero 
scalar. We therefore make all of the following analysis in 
T3. This means that we are making an analogous choice of 
tetrad in each metric, that all tetrad freedom is used up, and 
that ?/>04 represents a pure Type I effect. Indeed, it is clear 
from the definition 1331 that each of the ■!/'04 's in each of the 
three transverse symmetric tetrads is the square root of one 
of the three corresponding curvature invariants ('radiation 
scalars') introduced by Beetle & Burko (2002)^. However, 
only 7/)o4 in the quasi-Kinnersley frame vanishes in the 'Type 
D limit', and therefore its value is an invariant measure of 
the deviation from Type D. 

Finally, in any symmetric transverse frame (in partic- 
ular, in the quasi-Kinnersley tetrad T3) Eq. 12611 reduces 
to 



27[(- 



/ [(3^i 



2+^04)1 



(39) 



5* is equal to 1 in a Type D spacetime. (1 — 15*1) is also 
dimensionless; in this sense it is better than ■!/)04 as an in- 
variant measure of the deviation from Type D. We will refer 
to (1 — \S\) as the Non-Type-D-ness in the following. 



5.2 The Hartle-Thorne Weyl scalars in the 
quasi-Kinnersley frame 

The HT slow rotation parameter e is not a perturbative pa- 
rameter with respect to Type D: at linear order in e the HT 
solution in vacuum is equivalent to the Kerr metric, which of 
course is Type D at all orders in e. Discarding terms of 0{e^) 
and higher, the HT Weyl scalars in the quasi-Kinnersley 
(symmetric transverse, T3) tetrad have simple expressions: 



(V'o 



15 sin^ e 



it. 



X 

+ 

'2jHT ~ 



32MV3(r-2M)"^^~^^°"^ 
[3r^(r - 2Mf ln[r/{r - 2M)] 
2M{r - M){2M^ + 6rM - 3r^)] 



(40) 



M 3icos9J 
1 



16M6r7 

2 » ^4 



P2{ 



cos ( 



15r^ In 



- 2M 



lOMr^ (fiM'^ + 6M'\ + AM^r^ + 3Mr^ + 3r'^) 
16M'^(r - 6A/)] - 16MV(r - M)} 



+ 



1 



16A4"5r4 



P2(c 



15r* In 



2M 



- 10M(3r^ + 4rM^ + Sr'^M + 6M^)] Q . (41) 

For a Kerr BH Q — Qkctt = /M, and it is immediately 
seen that tpoi ~ (as it should). Writing the above relations 
as Taylor series of the form 



1p2 



(42) 



(2) 2 



+ 0(e 



it is easy to show that 



S = 3 



4'' 



' Tp2 



+ 0(e'). (43) 



Remarkably, the leading-order correction of S around 
its Type D value S = 1 is only of order e*: a fact due to 
S being a non-linear function of the Weyl scalars. Indeed, 
this doesn't mean that perturbations of order do not alter 
the speciality condition: rather, (1 — S) is quadratic in these 
perturbations*. This is quite a subtle point, that can better 
be understood by explicitly constructing the principal null 
directions. Such a construction has been worked out in detail 
by Cherubini et al. (2004) for the simpler case of the vacuum 
Kasner spacetimes: cf. in particular the discussion following 
their Eq. (2.15). Furthermore, at leading order in the ex- 
pansion the Non-Type-D indicators 1 — S and ipo4/il'2 are 
completely equivalent to each other. In Section|S|we show by 
explicit calculations that this equivalence of the two Type D 
indicators is in fact more general: it holds true (at least ap- 
proximately) for fully relativistic rotating star spacetimes as 
well. From Eq. 1401 we also see that {tpo4) ht is proportional 
to (Q — QKcrr). lu Other words, at leading order deviations 
from Type D are driven by the deviation of the quadrupole 



Although each of the three ipoi's is defined in a specific tetrad 
by Eq. 1331 . which is not per-se invariant under tetrad rotations 
of Type I and II, it can actually be expressed in terms of I and 
J only. See Beetle and Burko (2002) for more details. 



* It is easy to check, using an expansion like 1421 extended to 
higher order, that the leading fourth-order term in 1431 will not 
change. 
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moment from its Kerr value, as we would expect on physical 
grounds. 

From Eqs. H4()|l and l|41|l the asymptotic behavior of the 
Weyl scalars for large r at leading order is: 



(l/'04 
(V'2 



3sin^6'(Q-QKcr: 



2r5 



M 



(44) 
(45) 



Consequently, for the variables we use to measure deviations 
from Type D, at large r we have: 

3sin2(0) 



tp2 



s 



2Mr2 
27 sin'' (61) 



(Q — Qkc 



(0 



(46) 



(47) 



Finally, as a useful check of our results, we show that 
ip2 reduces to the correct limit in the Kerr case. The proof is 
quite straightforward. To avoid confusion, let us denote the 
standard Boyer-Lindquist coordinates for the Kerr metric as 
(s,a) instead of {r,9). Then the only non-zero Weyl scalar 
for the Kerr metric in the Kinnersley frame is given by 

— Af 

(^2) Kerr = 7 : Tg. (48) 

(s " la cos a)-^ 
The relation between the two sets of coordinates is^: 



a 

1 - 



1 + 
2M 



2M 



1 + 



1 - 
3M 



M 



' — a cos y sm t 



2r2 



r J 



(49) 



(50) 



Carrying out these coordinate transformations on equation 
14811 ■ setting J = Ma, Q — /M, and discarding terms of 
O(J^), we get 



(V'2)Ke 



M 



3i cos 9 J 



+ 2 



3 (r - M) 



(51) 



(6M + 5r) cos^ 9 - (2M + r)] / 



This same expression is obtained from equation I4H when 
Q = J'^/M is substituted. 

This confirms that in T3, with the parameters set for 
Kerr, we indeed have the Kinnersley frame, and that with 
these parameters, ip2 is the only non-zero Weyl scalar and 
matches the analytically known expression. Clearly the fur- 
ther reduction to Schwarzschild {J = Q — 0) also gives the 
correct known expression for (V'2)schw. 



5.3 The Manko and CST Weyl scalars in the 
quasi-KinnersIey frame 

For the Manko metric we computed the Weyl scalars in the 
quasi-Kinnersley frame by the procedure described in Sec- 
tion using the algebraic manipulation program Maple. 

^ For the HT metric to reduce to the Kerr metric we must set 
J = Ma, Q = J^/M and ignore terms of 0{J^). Formula (A5) in 
(Hartle &; Thorne 1968) contains a typo: in the r transformation, 
cos'^ 9 should be replaced by — cos^ 9. 



However the final analytic expressions are very long and 
unenlightening, so we do not give them here. As a check, 
we verified numerically that in the quasi-Kinnersley frame 
(V'04 )Manfco is zcro whcu the Manko parameters are set to 
Kerr (6^ = - M^) or Schwarzschild (a = 0, b'^ = -M^). 

For the CST metric we followed again the same pro- 
cedure, but in this case we had to solve a numerical prob- 
lem. The Weyl scalars depend on first and second order par- 
tial derivatives of the metric coefficients. For the analytical 
spacetimes (HT and Manko) the computation of the scalars 
is a trivial process. However, for the numerical spacetime one 
has to invoke numerical differentiation. The second order ra- 
dial derivatives turn out to be more sensitive to numerical 
error. Therefore we find that it is advantageous to avoid us- 
ing second order radial derivatives calculated straight from 
a second order finite difference scheme. Instead we use the 
Ricci identities in the Newman-Penrose formalism. These al- 
low us to express all second order radial derivatives in terms 
of second order angular derivatives and first order deriva- 
tives. Once we do this, the numerical calculation of the Weyl 
scalars turns out to be second order convergent on our nu- 
merical grid. 



6 NON-TYPE-D-NESS OF THE EXTERIOR 
VACUUM 

Using the procedure described in the previous Section we 
can compute the non-zero Weyl scalars tp2 and in the 
quasi-Kinnersley frame for all three metrics (HT, Manko and 
CST). Then we can evaluate the Non-Type-D-ness (1 - |S|) 
and the ratio |V'04/V'2|. To avoid unnecessary complications, 
in the following we will restrict consideration to the equato- 
rial plane. 

In Fig.l^we plot the invariant Non-Type-D-ness (1— jSj) 
as a function of the circumferential radius -Rcirc for the 
CST metric. Throughout this section we will always re- 
fer to circumferential radii, which are operationally well- 
defined quantities (whatever approximation we use for the 
metric in the stellar exterior). For example, for a HT model 
the circumferential stellar radius on the equator is simply 

r fHT) 1 

{Re, 9 = 7r/2) , where Re is the coordinate radius. 
The stellar radius in Manko coordinates pe can be evaluated 
using the procedure described in Section 4.1 of BS. Once 
again, it's a simple matter to compute the corresponding 

equatorial circumferential radius |^500''"''°' (Po ^ ~ 0)j for 
the Manko models. Of course the fact that we are using dif- 
ferent approximations for the metric will induce deviations 
between these stellar radii, but those deviations are typically 
below 1 % or so, even for the fastest rotating models. 

In both panels of Fig. lines from bottom to top give 
(1 — |S|) as a function of the circumferential radius -Rcirc 
for different values of the rotation rate. Since the calcula- 
tion only applies to the vacuum exterior spacetime, each 
line starts from the star's circumferential equatorial radius 
(Rcirc > Re)- The plot shows that (1 — \S\) has the essen- 
tial properties required of a Non-Type-D-ness: first of all, 
slower rotating model are supposed to be closer to Type D, 
and indeed (1 — 15*1) is smaller; secondly, as -Rcirc — > 00 the 
spacetime becomes asymptotically fiat and (1 — |S|) tends 
to zero. 
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Figure 5. Non-Type-D-ness (1 - |5|) for the FPS 14 (top panel) 
and FPS MM (bottom panel) sequences, computed using the CST 
metric, as a function of the circumferential radius -Rcirc- Within 
each panel, lines from bottom to top give (1 — |5|) as a function 
of the circumferential radius for different values of the rotation 
rate; slower rotating models are closer to Type D, so (1 — \S\) 
is smaller. Since the calculation only applies to the vacuum ex- 
terior spacetime, each line starts from the star's circumferential 
equatorial radius R^. 

For slowly rotating models (roughly speaking, for e < 
0.4), the lower lines in Fig.|^show that deviations from spe- 
ciality are at most of order 1 %. So in most astrophysical 
situations it may make sense to use a Type D approxima- 
tion: that is, to consider the exterior spacetime of a rotating 
star as a Kerr background, small corrections being induced 
by the star's oblate structure. This Type D approximation 
might be particularly useful for imposing boundary condi- 
tions in the computation of gravitational waves from rotat- 
ing relativistic stars. Indeed, the gravitational wave ampli- 
tude dies out only as 1/r, while our analysis shows that 
effects due to deviations from Type D die out faster than 
1/r. Therefore, far enough from the star the gravitational 
wave amplitude should (in some sense) be large enough that 
we can ignore the error made by using the Kerr approx- 



Figure 6. Close-up view of Fig. |21 Once again, the top panel 
refers to the FPS 14 sequence and the bottom panel to the FPS 
MM sequence. In red we overplot the Non-Type-D-ness for the 
Manko matching models, to visualize the difference with respect 
to the corresponding CST models. 



imation. On the other hand, the Kerr metric usually has 
a quadrupole moment that is 2 to 12 times smaller than 
that of a rapidly rotating star (Laarakkers & Poisson 1999). 
So, as long as we are in a regime in which the quadrupole 
moment affects the gravitational wave emission, the contri- 
bution of the quadrupole term will be underestimed by the 
corresponding factor. Only when we are far enough away 
that only the first-order terms in e matter do the neutron 
star and BH spacetimes agree. In this sense Kerr is not a 
good replacement for HT to order e'^, because of the large 
difference in the quadrupole moment. Using HT to order 
should give a much more accurate waveform consistent to 
this order. A wave extraction formalism which is consistent 
to second order in e is therefore worth developing. 

Fig. |S| shows a close-up view of the strong field region 
of Fig. |S] In red we overplot the Non-Type-D-ness we com- 
puted for the matching Manko models (when they exist: 
recall that the Manko solution can only be matched to fast 
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Figure 7. Ratio |V'04/i/'2l for the FPS 14 (top panel) and FPS 
MM (bottom panel) sequences, computed using the CST metric, 
as a function of the circumferential radius -Rcirc- See the caption 
of Fig.l^for more details. 



rotating models). For both sequences and for all values of 
the rotation rate the Manko metric is closer to Type D than 
the corresponding numerical model. 

In Fig.|7|we plot the ratio |i/'04/i/'2| as a function of the 
circumferential radius iicirc for the CST metric. As expected 
from our discussion in Sec. 15.11 the qualitative behavior is 
the same as in Fig. |^ \tp04/i^2\ is not invariant (^/>04 is a 
curvature invariant, in the sense discussed above, but i/>2 is 
not), but it has the advantage of being dimensionless. One 
may want to develop a 'quasi- Type-D' approximation for 
perturbations of rotating NSs in the framework of the NP 
formalism, similar to the Teukolsky formalism for BHs. In 
this context the perturbed ipo and i/>4 would represent ra- 
diation, typically decaying as 1/r. Then the value of the 
background tpo4 would be a useful number to consider for 
comparison. It is clear from Eq. 1441 for the HT approxi- 
mation and from Fig.|7|for the CST metric that ^04, decays 
faster than 1/r. 

One might think that the value of the Non-Type-D- 
ness at the equatorial stellar radius, which is a physically 



Figure 8. Non-Type-D-ness evaluated at the star's circumferen- 
tial equatorial radius i?e, as a function of Re, for the FPS 14 and 
FPS MM sequences (top and bottom panel, respectively). Filled 
circles correspond to CST models, empty circles correspond to 
matched HT models. 



well-identified location, could give an invariant measure of 
the deviation from speciality of each rotating stellar model. 
This is indeed the case for slowly rotating stars. In Fig. |H]we 
plot the value of (1 — 15*1) at the equatorial circumferential 
stellar radius i?e as a function of Re for our slowly rotat- 
ing FPS 14 and FPS MM sequences. This quantity has the 
behavior expected of a Non-Type-D-ness in this rotational 
range: both for the CST metric and for the HT metric it in- 
creases with Re (and with the rotation rate). Furthermore, 
it is larger for the FPS 14 sequence than for the FPS MM 
sequence. Once again, this is reasonable: low-mass models 
are more oblate than high-mass models, and their exterior 
spacetime is expected to deviate more from that of a BH. 

From Fig.l^and Fig.|S|it is clear that, at a given circum- 
ferential radius, (1 — |S|) increases with the stellar rotation 
rate. However, for large rotation rates (1 — \S\) at the stellar 
equator does not increase monotonically with rotation: we 
show this peculiar efi'ect in Fig. ^ In other words, this sin- 
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Figure 9. Non-Type-D-ness evaluated at the star's circumfer- 
ential equatorial radius _Re, as a function of iJe, for the two se- 
quences FPS 14 (circles) and FPS MM (squares). Filled symbols 
correspond to CST models, empty symbols correspond to match- 
ing Manko models (when they exist: cfr. BS). 



gle number cannot be used to measure the deviation from 
speciality of a very fast rotating stellar model. To under- 
stand the reason for this non-monotonic behavior we can 
look again at Fig. |S] At any given circumferential radius 
RciTc, (1 — IS'I) does increase as the star spins faster. How- 
ever the stellar circumferential radius (the value of i?circ at 
which each line starts) increases as well, due to the centrifu- 
gal deformation of the star. The two effects compensate, so 
that (1 — l^l) at the stellar radius shows a maximum as a 
function of the rotation rate. From Fig. |5]we see that this 
maximum shows up at the same rotation rate both for the 
CST metric and for the Manko matching models. 



7 CONCLUSIONS 

In this paper we compared three different approaches to the 
computation of rotating relativistic star spacetimes: the HT 
slow-rotation approximation, the Manko analytic vacuum 
solution, and numerical solutions of the full Einstein equa- 
tions obtained using the CST self-consistent field method. 
We first integrated the HT structure equations for five rep- 
resentative equations of state, keeping terms up to second 
order in the slow-rotation parameter e. Then we matched 
these models to the CST solutions, imposing that the grav- 
itational mass and angular momentum of the models be 
the same. We estimated the limits of validity of the slow- 
rotation expansion computing deviations in the spacetime's 
quadrupole moments and in the ISCO radii at different ro- 
tation rates. We found that deviations in the quadrupole 
moment are ~ 20 % for pulsars spinning with a period 
~ 1.5 milliseconds (the spin period of the fastest known pul- 
sar, PSR J1939-I-2134). However, for these same spin rates 
deviations in the ISCO radii are always smaller than ~ 1 %. 
Since the HT approximation gives excellent predictions for 
ISCOs up to the fastest pulsar spin periods, it can safely 



be used whenever a full numerical solution would be too 
cumbersome to implement. 

In the second part of the paper we focused on the ex- 
terior vacuum spacetime. We compared the HT approxima- 
tion and the Manko analytic solution with numerical mod- 
els using coordinate-independent quantities. For all three 
spacetimes we first introduced symmetric transverse frames, 
in which the only non-zero Weyl scalars are 7/12 and i/jq = 
ipi = '004 • From those frames we then selected what we call 
the 'quasi-Kinnersley' frame (cf. Nerozzi et al. 2004): the 
frame in which i/)2 is the only non-zero scalar when we set 
the parameters to reduce the metric to Type D. We com- 
puted 11)2 and '004 in this frame. The latter is a scalar cur- 
vature invariant (Beetle & Burko 2002) that vanishes for 
Type D, and therefore its value is a measure of what one 
would call 'Non-Type-D-ness'. For the HT metric the use 
of the quasi-Kinnersley frame led us to very simple analytic 
expressions for 7/12 and 004 • Due to their simplicity, our Eqs. 
1401 1 and 14111 may be useful in different contexts. Further- 
more, from Eq. 1401 we reached the important conclusion 
that, at leading order, the deviation from Type D is linear 
in the deviation of the quadrupole moment from its Kerr 
value, (Q — Qxerr). Therefore neutron stars with smaller 
mass and stiffer EOS deviate more from Type D. 

Finally we evaluated the (scalar curvature invariant) 
speciality index S introduced by Baker & Campanelli (2000) . 
Being dimensionless, (1 — |S|) is a better indicator of 'Non- 
Type-D-ness'. As i/)04, it increases with rotation rate and 
goes to zero at infinity for any given rotation rate. At leading 
order from Eq. 14311 (1 — S) oc {ipon / ip2f' ■ However, it is not 
easy to give a single number characterizing the Non-Type-D- 
ness of a given spacetime. One might think about using the 
value of (1 — 151) at the stellar radius for this purpose. This 
idea does not work in practice: since stellar models become 
more and more oblate as rotation increases, (1 — |S|) at the 
stellar radius eventually shows a maximum as a function of 
the rotation rate. 

Although deviations from Type D increase for fast rota- 
tion, our results show that they can be expected to be rather 
small (less than a few percent) for astrophysically relevant 
rotation rates. This suggests that, in perturbative calcula- 
tions of gravitational wave emission, rotating star exteriors 
can reasonably be approximated as Kerr-like spacetimes, or 
that one could develop an appropriate quasi- Type D ap- 
proximation. The tools we have illustrated here will help 
to quantify the accuracy of such approximations, and hence 
to estimate numerical errors involved in the calculations of 
gravitational waves from rotating relativistic stars. 

Acknowledgements 

We would like to thank Lior Burko for very kindly pro- 
viding us with notes on the method to identify the three 
transverse frames, and the anonymous referee for insight- 
ful comments. We also thank Marek Abramowicz, Joachim 
Almergren, Nils Andersson, Christian Cherubini, Wlodek 
Kluzniak, Kostas Kokkotas, John Miller, Andrea Nerozzi, 
Luciano RezzoUa, Adam Stavridis, Nick Stergioulas, Arun 
Thampan and Clifford Will for useful discussions. This 
work has been partly supported by the EU Programme 
'Improving the Human Research Potential and the Socio- 
Economic Knowledge Base' (Research Training Network 



Rotating neutron stars: an invariant comparison 15 



Contract HPRN-CT-2000-00137). FW is supported by a EP- 
SRC studentship. 



REFERENCES 

M. Abramowicz, G. J. E. Almergren, W. Kluzniak, A. V. Tham- 

pan, 2003, gr-qc/0312070 
N. Andersson, K. D. Kokkotas, 2001, Int. J. Mod. Phys. D 10, 

381. 

M. Ansorg, A. Kleinwacchter, R. Meinel, 2002, Astron. Astrophys. 
381, L49. 

D. C. Backer, S. R. Kulkarni, C. Heiles, M. M. Davis, W. M. Goss, 

1982, Nature 300, 615. 
J. Baker, M. CampancUi, 2000, Phys. Rev. D, 62, 127501. 

C. Beetle, L. M. Burko, 2002, Phys. Rev. Lett. 89, 271101. 

E. Berti, N. Stergioulas, 2004, MNRAS 350, 1416. 

M. Bocquet, S. Bonazzola, E. Gourgoulhon, J. Novak, 1995, As- 
tron. Astrophys. 301, 757. 

S. Bonazzola, E. Gourgoulhon, M. Salgado, J. -A. Marck, 1993, 
Astron. Astrophys. 278, 421. 

S. Bonazzola, E. Gourgoulhon, J. -A. Marck, 1998, Phys. Rev. D 
58, 104020. 

L. Burko, private communication, 2003. 

D. Chakrabarty, E. H. Morgan, M. P. Muno, D. K. Galloway, R. 
Wijnands, M. van der Klis, C. B. Markwardt, 2003, Nature 
424, 42. 

S. Chandrasekhar, 1983, The Mathematical Theory of Black 

Holes, Oxford University Press. 
S. Chandrasekhar, J. C. Miller, 1974, MNRAS 167, 63. 

C. Cherubini, D. Bini, M. Bruni and Z. Perjes, 2004, Class. Quant. 

Grav. 21, 4833. 

N. A. Collins, S. A. Hughes, 2004, Phys. Rev. D 69, 124022. 
G. B. Cook, S. L. Shapiro, S. A. Teukolsky 1994, Ap. J., 424, 
823. 

J. L. Friedman, J. R. Ipser, L. Parker, 1986, Ap. J. 304, 115. 
G. Fodor, Z. Perjes, 2000, Gen. Rel. Grav. 32, 2319-2332. 

G. Fodor, C. Hoenselaers, Z. Perjes, 1989, J. Math. Phys. 30, 
2252. 

J. B. Hartle, 1967, Ap. J., 150, 1005. 
J. B. Hartle, J. L. Friedman, 1975, Ap. J. 196, 653. 
J. B. Hartle, K. S. Thorne, 1968, Ap. J., 153, 807. 
K. D. Kokkotas, J. Ruoff, 2002, in Florence 2001, A relativistic 
spacetime odyssey, 325; see also gr-qc/0212105 

H. Komatsu, Y. Eriguchi, I. Hachisu, 1989, MNRAS, 237, 355. 
S. R. Kulkarni, 1995, in Millisecond Pulsars. A Decade of Sur- 
prise, volume 72 of ASP Conference Series, 79. 

W. G. Laarakkers, E. Poisson, 1999, Ap. J., 512, 282. 

V. S. Manko, E. W. Mielke, J. D. Sanabria-Gomez, 2000a, Phys. 

Rev. D 61, 081501(R). 
V. S. Manko, J. D. Sanabria-Gomez, O. V. Manko, 2000b, Phys. 

Rev. D 62, 044048. 
A. Nerozzi, C. Beetle, M. Bruni, L. M. Burko, D. PoUney, 

gr-qc/0407013 
J. C. Miller, 1977, MNRAS 179, 483. 

T. Nozawa, N. Stergioulas, E. Gourgoulhon, Y. Eriguchi, 1998, 
Astron. Astrophys. Suppl. Ser. 132, 431. 

D. PoUney, J. E. F. Skea, R. A. d'Inverno, 2000, Class. Quant. 
Grav. 17, 643. 

V. Re, M. Bruni, D. R. Matravers, F. T. White, 2003, Gen. Rel. 
Grav. 35, 1351. 

L. RezzoUa L., F. K. Lamb, S. L. Shapiro, 2000, Ap. J. 531, L141. 

F. D. Ryan, 1995, Phys. Rev. D 52, 5707. 
F. D. Ryan, 1997, Phys. Rev. D 55, 6081. 

M. Salgado, S. Bonazzola, E. Gourgoulhon, P. Haensel, 1994, As- 
tron. Astrophys. 291, 155. 
M. Shibata, M. Sasaki, 1998, Phys. Rev. D 58, 104011. 



N. R. SibgatuUin, 1991, Oscillations and waves in strong gravi- 
tational and electromagnetic fields (Engl, transl.), Springer, 
Berlin [orig. Russian, 1984, Nauka, Moscow]. 

H. C. Spruit, 1999, Astron. Astrophys. 341, LI. 

H. Stephani, 2004, Relativity: An Introduction to Special and 
General Relativity, Cambridge University Press. 

H. Stephani ,D. Kramer, M. MacCallum ,C. Hoenselaers ,E. Herlt, 
2003, Exact Solutions of Einstein's Field Equations, Cam- 
bridge University Press. 

N. Stergioulas, 2003, Rotating Stars in Relativity, Living Reviews 
in Relativity. 

N. Stergioulas, J. L. Friedman, 1995, Astrophys. J. 444, 306. 
J. M. Stewart, M. Walker, 1974, Proc. R. Soc. Lond. A 341, 49. 
M. Stute, M. Camenzind, 2002, MNRAS 336, 831. 
K. Sumiyoshi, J. Ma. Ibariez, J. V. Romero, 1999, Astron. Astro- 
phys. Suppl. Ser. 134, 39. 
S. Szekeres, 1965, J. Math. Phys. 6, 1387. 
S. A. Teukolsky 1973, Ap. J. 185, 635. 
F. Weber, N. K. Glendenning, 1992, Ap. J. 390, 541. 



16 Emanuele Berti, Frances White, Asimina Maniopoulou, Marco Bruni 

Table 1. Results of the HT integrations for EOS A. From left to 
right we list: the central energy density pc\ the stellar radius R, 
gravitational mass M and binding energy Eg for the non-rotating 
configuration; f2* = {M / R?)^/'^ . Then we give quantities of order 
e in rotation; the radius of gyration R* and the frame dragging 
function lo* at the center (normalized by H*). In the last five 
columns we give quantities of order e^, namely the rotational cor- 
rections on radius {SR* /R), mass (5M* /M) and binding energy 
{SEg/ Eb), the quadrupole moment Q* /{MR^) and the eccen- 
tricity e*. Units are given where appropriate. 



Pc (g cm ^) 


R (km) 


M/Mq 


Eg /Me 


Q* (km-i) 


r;/r 




~e2 
5R*/R 


SM*/M 


SE-^/Eb 


Q*/{MR^) 




1.25E15 


9.90 


1 


.04 


0.0933 


0.0398 


0, 


.583 


0.602 


0.120 


0.338 


0.417 


0.115 


1.13 


1.35E15 


9.86 


1 


.12 


0.111 


0.0416 


0, 


.592 


0.571 


0.114 


0.330 


0.395 


0.115 


1.12 


1.46E15 


9.81 


1 


.20 


0.129 


0.0434 


0, 


.601 


0.540 


0.107 


0.322 


0.372 


0.113 


1.11 


1.58E15 


9.74 


1 


.27 


0.148 


0.0451 


0, 


.609 


0.509 


0.100 


0.312 


0.347 


0.112 


1.11 


1.72E15 


9.66 


1 


.34 


0.168 


0.0468 


0, 


.617 


0.479 


0.0936 


0.302 


0.321 


0.110 


1.10 


1.86E15 


9.57 


1 


.40 


0.187 


0.0486 


0, 


.624 


0.450 


0.0870 


0.291 


0.295 


0.108 


1.08 


2.01E15 


9.47 


1 


.45 


0.205 


0.0503 


0, 


.631 


0.421 


0.0805 


0.280 


0.268 


0.107 


1.07 


2.18E15 


9.36 


1 


.50 


0.223 


0.0520 


0, 


.638 


0.394 


0.0742 


0.270 


0.242 


0.105 


1.06 


2.36E15 


9.25 


1 


.54 


0.238 


0.0537 


0, 


,644 


0.368 


0.0683 


0.260 


0.216 


0.103 


1.05 


2.56E15 


9.13 


1 


.58 


0.252 


0.0553 


0, 


.649 


0.344 


0.0628 


0.251 


0.191 


0.102 


1.04 


2.77E15 


9.00 


1 


.60 


0.264 


0.0570 


0, 


.654 


0.321 


0.0576 


0.242 


0.167 


0.101 


1.03 


3.00E15 


8.87 


1 


.62 


0.274 


0.0586 


0, 


.659 


0.300 


0.0528 


0.233 


0.144 


0.100 


1.02 


3.25E15 


8.75 


1 


.64 


0.281 


0.0601 


0, 


.663 


0.280 


0.0485 


0.225 


0.122 


0.0995 


1.01 


3.52E15 


8.62 


1 


.65 


0.286 


0.0617 


0, 


.667 


0.261 


0.0445 


0.218 


0.100 


0.0988 


1.00 


3.81E15 


8.49 


1 


.65 


0.289 


0.0632 


0, 


.670 


0.244 


0.0408 


0.211 


0.0797 


0.0982 


0.995 


4.13E15 


8.36 


1 


.66 


0.290 


0.0647 


0, 


.673 


0.228 


0.0376 


0.205 


0.0601 


0.0977 


0.987 



Table 2. Same as in TablelT] but for EOS AU. 



Pc (g cm ^) 


R (km) 


M/Mq 


Eb/Mq 


Q* (km-i) 


r;/r 




SR*/R 


SM*/M 


SE-^/Eb 


Q*/{MR^) 




7.62E14 


10.3 


0.602 


0.0254 


0.0285 


0.513 


0.763 


0.150 


0.322 


0.461 


0.0946 


1.12 


8.36E14 


10.3 


0.742 


0.0413 


0.0317 


0.539 


0.718 


0.142 


0.342 


0.468 


0.107 


1.13 


9.16E14 


10.3 


0.895 


0.0632 


0.0347 


0.562 


0.670 


0.134 


0.351 


0.460 


0.115 


1.14 


1.00E15 


10.3 


1.05 


0.0914 


0.0375 


0.582 


0.620 


0.124 


0.350 


0.439 


0.119 


1.14 


1.10E15 


10.4 


1.22 


0.128 


0.0402 


0.601 


0.566 


0.112 


0.342 


0.406 


0.120 


1.13 


1.21E15 


10.4 


1.40 


0.175 


0.0429 


0.621 


0.508 


0.0989 


0.327 


0.361 


0.120 


1.12 


1.32E15 


10.4 


1.57 


0.229 


0.0456 


0.639 


0.450 


0.0850 


0.310 


0.311 


0.119 


1.10 


1.45E15 


10.3 


1.72 


0.285 


0.0480 


0.656 


0.397 


0.0716 


0.291 


0.260 


0.118 


1.08 


1.59E15 


10.3 


1.84 


0.338 


0.0502 


0.671 


0.349 


0.0594 


0.273 


0.210 


0.118 


1.06 


1.74E15 


10.2 


1.94 


0.383 


0.0522 


0.684 


0.309 


0.0489 


0.258 


0.165 


0.118 


1.04 


1.91E15 


10.1 


2.01 


0.421 


0.0540 


0.694 


0.274 


0.0399 


0.244 


0.124 


0.118 


1.03 


2.09E15 


9.94 


2.06 


0.451 


0.0556 


0.703 


0.244 


0.0323 


0.232 


0.0873 


0.119 


1.01 


2.29E15 


9.81 


2.10 


0.473 


0.0573 


0.711 


0.217 


0.0256 


0.222 


0.0533 


0.120 


1.00 


2.51E15 


9.68 


2.12 


0.488 


0.0588 


0.717 


0.194 


0.0199 


0.213 


0.0226 


0.122 


0.990 


2.76E15 


9.54 


2.13 


0.497 


0.0602 


0.722 


0.174 


0.0150 


0.205 


-0.00527 


0.123 


0.980 


3.02E15 


9.40 


2.13 


0.499 


0.0616 


0.727 


0.156 


0.0110 


0.198 


-0.0305 


0.123 


0.971 
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Table 3. Same as in Table[I] but for EOS FPS. 



Pc (g cm 3) 


R (km) 


M/Mq 


Eb/Mq 


Q* (km-i) 


r;/r 




SR'/R 


SM'/M 


5E*g/EB 


Q*/{MR2) 


e* 


8.Uvitjl4 


11.2 


0.873 


0.0533 


0.0303 


0, 


.546 


0.695 


0.139 


0.339 


0.454 


0.108 


1.13 


8.84E14 


11.1 


0.986 


0.0704 


0.0324 


0, 


.560 


0.659 


0.132 


0.339 


0.441 


0.111 


1.13 


9.73E14 


11.1 


1.10 


0.0901 


0.0345 


0, 


.573 


0.621 


0.125 


0.335 


0.422 


0.112 


1.13 


1.07E15 


11.0 


1.21 


0.112 


0.0364 


0, 


.585 


0.584 


0.117 


0.328 


0.399 


0.112 


1.12 


1.18E15 


10.9 


1.31 


0.135 


0.0384 


0, 


.595 


0.547 


0.109 


0.318 


0.372 


0.111 


1.11 


1.30E15 


10.8 


1.40 


0.159 


0.0403 


0, 


.605 


0.512 


0.101 


0.308 


0.344 


0.109 


1.10 


1.43E15 


10.7 


1.48 


0.183 


0.0421 


0, 


.614 


0.477 


0.0935 


0.296 


0.315 


0.107 


1.09 


1.57E15 


10.6 


1.56 


0.206 


0.0439 


0, 


.622 


0.444 


0.0861 


0.285 


0.286 


0.105 


1.08 


1.73E15 


10.5 


1.62 


0.227 


0.0457 


0, 


.630 


0.412 


0.0790 


0.273 


0.256 


0.103 


1.07 


1.91E15 


10.3 


1.67 


0.247 


0.0475 


0, 


.637 


0.382 


0.0721 


0.261 


0.227 


0.102 


1.05 


2.10E15 


10.1 


1.71 


0.264 


0.0493 


0, 


.643 


0.353 


0.0657 


0.250 


0.198 


0.0999 


1.04 


2.31E15 


9.97 


1.75 


0.278 


0.0510 


0, 


.649 


0.327 


0.0598 


0.240 


0.170 


0.0985 


1.03 


2.54E15 


9.80 


1.77 


0.289 


0.0528 


0, 


.654 


0.302 


0.0543 


0.230 


0.144 


0.0974 


1.02 


2.80E15 


9.62 


1.79 


0.297 


0.0545 


0, 


.658 


0.279 


0.0495 


0.221 


0.118 


0.0964 


1.01 


3.08E15 


9.45 


1.80 


0.301 


0.0561 


0, 


.662 


0.258 


0.0451 


0.213 


0.0941 


0.0955 


0.999 


3.39E15 


9.28 


1.80 


0.303 


0.0577 


0, 


,666 


0.239 


0.0413 


0.206 


0.0713 


0.0947 


0.990 



Table 4. Same as in Table but for EOS L. 













~ e 




~e2 










Pc (g cm ^) 


R (km) 


M/Mq 


Eb/Mq 


Q* (km-i) 


R*/R 


uj*/n' 


5R*/R 


SM*/M 


5E%/Eb 


Q*/{MR^) 


e* 


2.35E14 


14.9 


0.362 


0.00513 


0.0127 


0.426 


0.886 


0.173 


0.234 


0.396 


0.0522 


1.07 


2.65E14 


14.5 


0.502 


0.0119 


0.0156 


0.473 


0.853 


0.164 


0.302 


0.459 


0.0804 


1.11 


3.00E14 


14.5 


0.680 


0.0242 


0.0182 


0.510 


0.812 


0.159 


0.350 


0.500 


0.105 


1.14 


3.39E14 


14.6 


0.896 


0.0445 


0.0206 


0.540 


0.765 


0.153 


0.377 


0.517 


0.122 


1.16 


3.83E14 


14.8 


1.14 


0.0751 


0.0228 


0.564 


0.711 


0.144 


0.387 


0.510 


0.132 


1.17 


4.33E14 


15.0 


1.40 


0.116 


0.0249 


0.585 


0.655 


0.134 


0.383 


0.485 


0.135 


1.16 


4.89E14 


15.1 


1.67 


0.168 


0.0268 


0.603 


0.597 


0.121 


0.370 


0.448 


0.133 


1.15 


5.53E14 


15.1 


1.91 


0.229 


0.0286 


0.619 


0.540 


0.108 


0.351 


0.402 


0.130 


1.14 


6.24E14 


15.1 


2.14 


0.292 


0.0303 


0.633 


0.486 


0.0955 


0.330 


0.353 


0.126 


1.12 


7.06E14 


15.0 


2.32 


0.353 


0.0319 


0.645 


0.437 


0.0834 


0.309 


0.304 


0.122 


1.10 


7.97E14 


14.8 


2.46 


0.407 


0.0334 


0.656 


0.393 


0.0727 


0.290 


0.259 


0.118 


1.08 


9.01E14 


14.6 


2.56 


0.450 


0.0348 


0.664 


0.355 


0.0636 


0.273 


0.218 


0.116 


1.06 


1.02E15 


14.4 


2.64 


0.483 


0.0361 


0.671 


0.322 


0.0560 


0.258 


0.182 


0.113 


1.05 


1.15E15 


14.2 


2.68 


0.504 


0.0373 


0.676 


0.294 


0.0498 


0.246 


0.151 


0.111 


1.03 


1.30E15 


13.9 


2.70 


0.514 


0.0383 


0.679 


0.271 


0.0449 


0.235 


0.125 


0.110 


1.02 


1.47E15 


13.7 


2.70 


0.517 


0.0393 


0.682 


0.251 


0.0411 


0.226 


0.102 


0.108 


1.01 
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Table 5. Same as in Table[I] but for EOS APRb. 



Pc (g cm 3) 


R (km) 


M/Mq 


Eb/Mq 


Q* (km-i) 


r;/r 




SR*/R 


SM'/M 


5E*g/EB 


Q*/{MR2) 


e* 


6.16E14 


12.1 


0.654 


0.0243 


0.0234 


0.497 


0.771 


0.152 


0.294 


0.430 


0.0835 


1.10 


6.78E14 


11.9 


0.776 


0.0368 


0.0261 


0.518 


0.733 


0.145 


0.311 


0.438 


0.0927 


1.11 


7.46E14 


11.8 


0.917 


0.0546 


0.0288 


0.539 


0.691 


0.137 


0.322 


0.436 


0.100 


1.12 


8.21E14 


11.7 


1.07 


0.0787 


0.0315 


0.559 


0.644 


0.128 


0.326 


0.423 


0.105 


1.12 


9.04E14 


11.6 


1.24 


0.110 


0.0341 


0.578 


0.594 


0.118 


0.324 


0.400 


0.108 


1.12 


9.95E14 


11.6 


1.40 


0.147 


0.0366 


0.596 


0.543 


0.107 


0.316 


0.369 


0.109 


1.11 


1.10E15 


11.5 


1.56 


0.190 


0.0390 


0.613 


0.492 


0.0953 


0.305 


0.331 


0.109 


1.10 


1.21E15 


11.4 


1.70 


0.236 


0.0413 


0.629 


0.443 


0.0837 


0.291 


0.289 


0.109 


1.08 


1.33E15 


11.3 


1.83 


0.282 


0.0435 


0.643 


0.397 


0.0724 


0.276 


0.247 


0.108 


1.07 


1.46E15 


11.1 


1.94 


0.325 


0.0456 


0.656 


0.355 


0.0619 


0.262 


0.204 


0.108 


1.05 


1.61E15 


11.0 


2.02 


0.365 


0.0475 


0.668 


0.316 


0.0522 


0.248 


0.162 


0.108 


1.03 


1.77E15 


10.8 


2.09 


0.398 


0.0494 


0.678 


0.281 


0.0435 


0.235 


0.123 


0.108 


1.02 


1.95E15 


10.6 


2.14 


0.425 


0.0512 


0.688 


0.250 


0.0359 


0.224 


0.0865 


0.109 


1.01 


2.15E15 


10.5 


2.17 


0.444 


0.0529 


0.695 


0.223 


0.0292 


0.214 


0.0525 


0.110 


0.993 


2.36E15 


10.3 


2.19 


0.457 


0.0545 


0.702 


0.199 


0.0235 


0.205 


0.0212 


0.111 


0.981 


2.60E15 


10.1 


2.20 


0.463 


0.0560 


0.707 


0.178 


0.0187 


0.198 


-0.00756 


0.112 


0.971 
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Table 6. Properties of HT models matching the slowest rotating 
CST models listed in Tables 1-5 of BS. For each EOS, on top we 
present matching models for the M = 1.4M© ('14') sequence; 
separated by an horizontal line, we give matching models for the 

maximum-mass ('MM') sequence. From left to right we list: the 
central energy density pc for which we find a matching solution 
(if such a central energy density exists); the value of e = J/ J* at 
the given central energy density; the gravitational mass M and 
the angular nioiiientum J, in geometrical units; the value of the 
quadrupole moment predicted by the HT equations of structure, 
Q^"^ , in geometrical units; the percentage deviation of the HT 
quadrupole moment from its value Q for the CST spacetime. 



EOS 


n„ ( p" cvr\ \ 
He \s ^^^^ ) 




M (km) 


J (km-^) 




SO ('%') 


A 14 


1.8110E15 


0.23236 


2.0729 


0.8121 


1.1053 


10.4 




1.7238E15 


0.40060 


2.0800 


1.3267 


3.2644 


22.9 




1.6152E15 


0.55642 


2.0870 


1.7037 


6.2084 


41.8 




J.. 4141 1 OCjIO 


U. 1 \JOaO 




9 n99n 


1 1 •^TO 
ll.O 1 z 


HA 9 


A MM 


3.4253E15 


0.19954 


2.4519 


0.9857 


0.7192 


15.4 




3.0369E15 


0.32653 


2.4608 


1.5731 


2.0107 


21.3 




2.6524E15 


0.46715 


2.4726 


2.1471 


4.2850 


32.4 




2.2177E15 


0.65300 


2.4870 


2.7173 


8.6795 


59.7 


AU 14 


1.1916E15 


0.24345 


2.0734 


0.8704 


1.5633 


10.4 




1.1666E15 


0.40245 


2.0794 


1.3551 


4.1490 


23.3 




1 191 ypi ^ 


u. oyyyo 


Z.UoDo 




Q fiQI 7 

o.oyi 1 






1.0426E15 


0.87088 


2.0925 


2.0808 


16.160 


109.8 


AU MM 


2.5557E15 


0.21819 


3.1645 


1.9389 


1.6913 


23.5 




2.2777E15 


0.36512 


3.1838 


3.1382 


4.7843 


29.3 




2.0247E15 


0.51420 


3.2068 


4.1755 


9.4664 


40.5 




1.7585E15 


0.69420 


3.2302 


5.0587 


16.844 


66.1 


FPS 14 


1.2642E15 


0.25144 


2.0720 


0.8722 


1.6662 


10.2 




1.2195E15 


0.39258 


2.0769 


1.2971 


4.0094 


21.1 




1.1565E15 


0.54440 


2.0825 


1.6649 


7.5270 


39.0 




1.0626E15 


0.74208 


2.0881 


1.9766 


13.279 


74.8 


FPS MM 


2.8003E15 


0.20109 


2.6648 


1.1604 


0.9524 


14.7 




2.4765E15 


0.32814 


2.6740 


1.8467 


2.6577 


20.8 




2.1574E15 


0.46924 


2.6864 


2.5193 


5.6832 


32.1 




1.7968E15 


0.65491 


2.7014 


3.1853 


11.536 


59.7 


L 14 


4.2540E14 


0.25862 


2.0717 


0.9726 


4.0544 


10.4 




4.1435E14 


0.43170 


2.0760 


1.5007 


10.756 


27.1 




3.9803E14 


0.61954 


2.0802 


1.8991 


20.412 


54.6 


L MM 


1.2480E15 


0.18716 


4.0130 


2.5589 


3.0214 


13.8 




1.0717E15 


0.32762 


4.0297 


4.3607 


9.6976 


21.0 




9.2990E14 


0.47634 


4.0515 


6.0142 


21.145 


34.0 




7.8669E14 


0.66069 


4.0754 


7.4955 


41.295 


63.6 


APRb 14 


9.7905E14 


0.27032 


2.0753 


0.9349 


2.1641 


10.1 




9.6661E14 


0.35913 


2.0781 


1.2034 


3.7615 


17.4 




9.5231E14 


0.44237 


2.0810 


1.4276 


5.6017 


25.5 




9.3547E14 


0.52688 


2.0838 


1.6238 


7.7629 


35.6 




9.1459E14 


0.62069 


2.0866 


1.8017 


10.442 


49.1 


APRb MM 


2.4553E15 


0.17324 


3.2640 


1.6051 


1.1343 


25.3 




2.3001E15 


0.24355 


3.2705 


2.2296 


2.2707 


23.0 




2.0765E15 


0.36617 


3.2864 


3.2541 


5.2083 


28.3 




1.8261E15 


0.52120 


3.3086 


4.3615 


10.632 


39.8 




1.6702E15 


0.63017 


3.3232 


4.9670 


15.494 


54.0 




1.4393E15 


0.81746 


3.3393 


5.5922 


25.448 


94.7 



